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A COMPILATION OF CURRENT COMPUTER PROGRAMS 


FOR LOW-THRUST TRAJECTORY ANALYSIS 

By Alfred C. Mascy 

NASA Headquarters 
Moffett Field, Calif. 

SUMMARY 


This paper provides a catalogue of brief descriptions of current low- 
thrust trajectory and mass computation programs. All known sources were 
contacted, and questionnaires were sent to those who indicated that they pos- 
sessed such codes. A total of 59 descriptions of computer programs pertinent 
to this discipline are compiled. They are capable of generating flyby and 
rendezvous trajectories with their associated mass computation. Included in 
the codes are those capable of variable thrust trajectory as well as constant 
thrust with coast trajectory calculations. No attempt is made to evaluate 
any computer code or in any way compare programs; rather, the descriptions of 
existing codes are presented. Further details concerning a particular pro- 
gram may be obtained from the organization shown on each description. Many 
programs are well documented, and all organizations are encouraged to use the 
services provided by the regional dissemination center, COSMIC, located at 
the University of Georgia and sponsored by the NASA Office of Technology 
Utilization. 


INTRODUCTION 


There is increasing interest in the application of advanced propulsion 
systems, such as electric rockets, to advanced missions having high velocity 
requirements. Better utilization of existing computer programs, or the devel 
opment of new codes supplemented by experience with existing programs, is 
needed to define these requirements. This survey indicates that numerous 
low-thrust computer programs now exist, but many are not reported in the 
literature. The majority of these programs are worthwhile and can be of 
greater use if they are made known to the community of low-thrust analysts. 

Thus, the purpose of this paper is to provide a catalogue of brief 
descriptions of current low-thrust trajectory computer programs. This paper 
does not attempt to evaluate any computer code or to compare one program with 
another. The descriptions of computer codes pertinent to low-thrust mission 
analysis serve only to introduce the respective programs in order to advance 
research in this area by the exchange of ideas, trajectory concepts, and 
mathematical methods already applied. 

Gratitude is expressed to all organizations participating in this survey 
they gave freely of their time and information to enhance the scope of this 
work . 



SYMBOLS 


a Q initial acceleration 

Isp specific impulse of thruster 

T 

thrust-to-weight ratio 
Voo hyperbolic excess speed 

a power-plant specific mass 

n efficiency of thrusters 

payload mass fraction 
Vtp propellant mass fraction 

y w power-plant mass fraction 

COMPILATION OF PROGRAM DESCRIPTIONS 


All organizations that might possess a low-thrust computer program 
capability were contacted. After an initial contact by telephone, question- 
naires (see sample on page 7) were sent to the individuals who wished to par- 
ticipate in the survey. The data on each completed questionnaire were 
transferred to an information sheet that did not include the multiple choice 
or nonapplicable answers. Some questionnaires were only partially completed 
but were accompanied by the program documentation, which was transferred to 
the information sheet. It is hoped that this transfer was done as accurately 
as the programmer would have done it. Copies of the information were sent to 
all participants for their comments, corrections, and updating. The present 
publication reflects all changes. 

A total of 59 descriptions of programs pertinent to low- thrust 
trajectory and mass computation are compiled (pp. 8-66). It is believed that 
these represent most of the codes in this discipline. These descriptions 
represent approximately 52 different programs; many have been modified and are 
described as slightly different codes by various users. Rather than screen a 
program that may have been modified for the user's needs, all low-thrust code ' 
descriptions submitted have been included in this compilation. However, some 
programs were not directly applicable to low-thrust trajectory analysis and 
have not been included in this survey. 

There are about twice as many two-dimensional codes as there are three- 
dimensional codes. Most two-dimensional codes utilize a polar coordinate 
system, while most three-dimensional programs use a Cartesian system. Only 
two programs actually appear to account for more than the usual two bodies in 
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a gravitational system. Most codes have been programmed in Fortran IV, 
although a few require a knowledge of machine language. Most codes do require 
initial guesses such as the Lagrange multipliers; however, a few programs 
generate the necessary starting values from an approximate solution, such as 
a circular arc trajectory or an iterated circle start solution. A few pro- 
grams are applicable only to the planetocentric phase, and the others are 
flyby and rendezvous decks; most programs are capable of both operations. 

There is an almost total absence of round-trip trajectory codes. Approxi- 
mately 15 percent of the codes include some estimate of the planetocentric 
phase, as well as the heliocentric phase. 

Potential users who are interested in further information on the programs 
are advised to contact the respective company. Many programs are well docu- 
mented, and all organizations are encouraged to submit them to COSMIC for 
inclusion in the library of available programs. COSMIC is one of the newest 
regional dissemination centers located at the University of Georgia and spon- 
sored by the NASA Office of Technology Utilization. In general, its purpose 
is to acquire and disseminate information and software (card decks, tapes, 
etc.) on technical and scientific computer programs originated during NASA 
work. Once a program is accepted by COSMIC, the authoring organization no 
longer needs to prepare and transmit its program to potential users. Organi- 
zations are encouraged to use the dissemination service. 
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Sample Questionnaire 

Aerojet General, Space 
General Plant 

AFIT-SEE, Wright-Patterson 
AFB 

Analytical Mechanics Associ- 
ates 

Analytical Mechanics Associ- 
ates 

AVCO 

Bell Telephone Laboratories 

Boeing/Space Division 

Boeing/Space Division 

Colorado State University 

Computer Applications, Inc. 

Connecticut University, Aero- 
space Engineering Dept. 

General Dynamics /Convair 

General Dynamics/Fort Worth 

General Dynamics/Fort Worth 

General Dynamics/Fort Worth 

General Dynamics Fort/Worth 

General Electric Company 

General Electric - M§SD 

Goddard Space Flight Center, 
NASA 

Goddard Space Flight Center, 
NASA 

Grumman Aircraft 


Name of Program 


ORPERT 

Low-Thrust Solar Probe 

2-Body Trajectory Optimization 
Program 

Second Variation Rocket Trajec- 
tory Optimization Program 

#1546 - Generalized Variation 
of Parameters Program 

DEPAK; MULTI 

Variable Thrust Polynomial 
Optimizer 

Space Trajectory Optimization 
Program (AS 2080) 

Characteristic Length Method 

PDER, PDERDP (double precision 
version) 

Convair Low-Thrust Optimization 
Programs 

Y28 Optimum Low-Thrust Procedure 

Low-Thrust Trajectory Optimization 

Interplanetary Low-Acceleration 
Constant-Thrust Trajectory 
Patching Program 

Low-Thrust Propulsion System 
Optimization Program 

POLAR 

Finite Burn Orbit Insertion 
Program 

Variable Power, Variable 
Thrust 0ption/F369 

Low-Thrust Trajectory Subprogram 
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Low-Thrust Interplanetary 
Trajectory Optimization 



Page 

Comp any/Organization 

Name of Program 

29 

Hayes International Corp. 

Orbital Elements for Low-Thrust 
Trajectories 

30 

Hayes International Corp. 

Quasi-Steady Orbital Rendezvous 

31 

Hughes Aircraft Company 

Low-Thrust Trajectory 

32 

IBM Federal Systems Division 

OPGUIDE 

33 

IIT Research Institute 

LTNAV-1 (Low-Thrust Navigation 
Analysis) 

34 

I IT Research Institute 

LTNAV-2 (Low-Thrust Planetocentrii 
Navigation Analysis) 

35 

JPL 

Optimum Thrust Programs for Power 
Limited Propulsion Systems 
(1016000) 

36 

JPL 

ASTRAL 

37 

JPL 

MALTS 

38 

Langley Research 
NASA 

PRESTO (DELTA) 

39 

Langley Research Center, 
NASA 

TOMSD (Trajectory Optimization 
by Method of Steepest Descent) 

40 

Lewis Research Center, NASA 

Lewis N-Body Code 

41 

Lewis Research Center, NASA 

MAXIM 

42 

Lockheed Missiles and Space 
Company, Research Labora- 
tories 

Low-Thrust Heliocentric 
Trans fers 

43 

Lockheed Missiles § Space 
Company, Research Labora- 
tories 

Low -Thrust Escape From Planet 

44 

Lockheed Missiles § Space 
Company 

LTC0MP 

45 

Los Alamos Scientific 
Laboratory 

1965 Version, Produced by Dept, 
of Aerospace and Mechanical 
Sciences, Princeton University, 
Report No. 7170-1 

46 

McDonnell Douglas 
Astronautics 

FYNDIF 

47 

McDonnell Douglas 
Astronautics 

LOW TOP 

48 

McDonnell Douglas 
Astronautics 

WHAMO 

49 

McDonnell Douglas 

HIT0P 


Astronautics 
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Page 

Company/Organi zation 

Name of Program 

50 

North American Rockwell 
Corp., Space Division 

Optimum Trajectory Program for 
Power-Limited Propulsion 
Systems 

51 

OART, Mission Analysis 
Division, NASA 

BABE 

52 

OART, Mission Analysis 
Division, NASA 

QUICKLY 

53 

Princeton University, ASMAR 

Lion 1 

54 

Princeton University, ASMAR 

Campbell 1 (modified item) 

55 

Princeton University, ASMAR 

Gordon 2 

56 

Princeton University, ASMAR 

PRIMER 

57 

Rand Corporation 

ROCKET 

58 

RCA/ Astro Electronics 
Division 

Gordon 2 

59 

System Sciences Corp. (sub- 
sidiary of Computer 
Science Corp.) 

High-Precision, Low-Thrust 
Integration and Optimization 
Program 

60 

Texas University 

NMS 

61 

TRW Systems 

MALTS, Mission Analysis & 
Simulation Program 

62 

United Aircraft Research 
Laboratories 

Mass Ratio Optimization (sim- 
plified and improved payload 
fraction definition) F487 

63 

United Aircraft Research 
Laboratories 

High and Low Thrust Mass 
0ptimization/F211 

64 

United Aircraft Research 
Laboratories 

Mass Ratio Optimization (sim- 
plified payload fraction) F530 

65 

United Aircraft Research 
Laboratories 

Variable Power, Variable Thrust 
0ption/F369 

66 

United Aircraft Research 
Laboratories 

Constant -Thrust, Multiple- 
Coast/F615 
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LOW- THRUST TRAJECTORY AND MASS COMPUTATION PROGRAM 


Company Organization 
Name of Program 
Authors/Contact 


Program Funded by 


i.e., company, NASA contract, . . . 
Description 


Number of Dimensions 
Number of Bodies 


Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 

Other Options 
Input (briefly) 


(Please circle or add appropriate comment.) 
sions : i.e. , 1, 2, 3 

s : i.e. , 2, 3, . . . 

em : i.e., Cartesian, polar, . . . 

bnique : i.e., variable step, fixed step. 


i.e., Fortran IV, II, . . . 
i.e., IBM 7094, CDC 6600, Univac 1108, . . . 
i.e., flyby, rendezvous, round trip, . . . 
i.e., variable, constant, . . . 
i.e., constant, function of time, function of 
distance, .. . . 

i.e., position and velocity,' Julian dates, hyper- 
bolic excess velocity, ISP, powerplant specific 


mass, a, . . . 

Initial Guesses Required: i.e., Lagrange multipliers, T/W, . . . 

Parameters Optimized : i.e., ISP, T/W 0 , V OT , power level, /a 2 dt, 

coast time, thrust angle, . . . 

Optimization Technique : i.e., gradient method, Newton-Raphson, . . . 

Output (briefly) : i.e., time history, ISP, T/W, payload, /a 2 dt,. . . 

Planetocentric Phase : not included in program, high-thrust capture and 

departure, low-thrust spiralling, . . . 

Limitations of Program : i.e., travel angle, T/local g, . . . 

Other Applications of Program 

Documentation of Program 

Availability to Other Users 

Other Comments 

Please Mail Completed Form to: Alfred C. Mascy 

NASA Mission Analysis Division, OART 
Mail Stop 202-8 

Moffett field, California 94035 
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Company/Organization 

Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 

Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Input (briefly) 


Output (briefly) 
Limitations of Program 
Other Applications 


Documentation 
Availability 
Other Comments 


Aerojet General, Space General Plant, 

El Monte, California 
: ORPERT 

: Dr. Leonard Pode/Wayne F. Brady 

: Company 

Description 

: 3 

: Up to 5 satellites and Earth, Sun and Moon 

currently 
: Cartesian 

: Fixed step 

: Fortran IV 

: IBM 360/40 

: Currently written for Earth orbit 

: Constant 

: Position and velocity, Julian dates, time to 

start § stop thrust, direction of thrust WRT 
VEL vector, magnitude of thrust, weight flow 
rate 

: Time history 

: Eccentricity / 0 

: Will handle up to 5 satellites in Earth orbit; 

includes perturbations due to atmosphere. 
Earth harmonics. Moon § Sun 
: A report is in rough draft. 

: Program not completed due to funding. 

JPL ephemeris tape used. Conversion to another 
gravitational center would not be difficult. 
Program originally written to simulate 
deployment and separation history of several 
satellites in orbit. 
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Company/Organization 
Name of Program 
Authors /Contact 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Limitations of Program 
Other Applications 
Documentation 
Availability 


AFIT-SEE, W-Patterson AFB, Ohio 45433 
Low -Thrust Solar Probe 

Major Roger W. Johnson/Capt . H. T. Brock 
Description 

2 

2 

Polar 

Variable step 
Fortran IV 
IBM 7094 
Flyby 
Constant 

Position and velocity 

T/W 

/a 2 dt 

Newton-Raphson 
Time history 
Not included in program 
Travel angle 

Available approximately June 1968 
Thesis form 
Yes (listing) 



C omp any / Or g an i z at i on 
Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 


Initial Guesses Required 

Parameters Optimized 
Optimization Technique 


Planetocentric Phase 

Documentation 

Availability 


Analytical Mechanics Associates, Inc. 

2-Body Trajectory Optimization Program 

J. H. Campbell 

Sub -contract from Princeton 

Description 

3 

N 

Cartesian 
In mass ratio 
Fortran IV 
IBM 360 

Flyby, rendezvous 
Constant 

Constant, function of distance 
Position and velocity, Julian dates, hyperbolic 
excess velocity, ISP, powerplant specific 
mass, alpha, launch vehicle, sizing 
parameters, etc. 

Lagrange multipliers, T/W, and any unspecified 
initial conditions 

ISP, T/W 0 , \f m , coast time, thrust angle, t G , tf 
Employs indirect method coupled with gradient 
method for solution of boundary value 
problem 

High -thrust departure 
In progress 

Available at Princeton U. and GSFC 
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Company/ Organization 
Name of Program 

Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Limitations of Program 

Documentation 

Availability 

Reference 


Analytical Mechanics Associates, Inc. 

Second Variation Rocket Trajectory 
Optimization Program 

Dr. Henry J. Kelley, Westbury, N. Y. Office 
Various NASA MSC contracts 

Description 

3 

2 

Cartesian 

Fixed step Runge-Kutta 
Fortran IV 
IBM 360 

Rendezvous and orbit transfer options 
Constant 

Constant or solar 
Initial and terminal conditions 
First guess on steering history 
Steering histories and switch times 
Second variation 

Time history of control, state and multipliers 
Not included 

Program is primarily hi ght- thrust, but should 
hqndle low- thrust heliocentric well. 

Presently undocumented 
Has been used exclusively in-house 
"A Trajectory Optimization Technique Based Upon 
the Theory of the Second Variation," 

H. J. Kelly, R. E. Kopp and H. G. Moyer, 

AIAA No. 63-415, Aug. 19, 1963. 


11 



Company/Organization 
Name of Program 

Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 

Program Language 
Computer Now Used 
Trajectory Options 

Specific Impulse 
Power Options 
Other Options 


Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Output (briefly) 
Planetocentric Phase 

Other Applications 
Documentation 
Availability 
Other Comments 


AVCO 

#1546 - Generalized Variation of Parameters 
Program 

F. H. Scammell, M. Curley/D. P. Fields 
Company 

Description 

3 

N-Bodies (all planets, plus Sun and Moon) 
Cartesian 

Variable step (Adams 4th order predictor- 
corrector) 

Fortran IV 

IBM 360/75, 360/65 

Central body determination and switchover, 
flyby, round trip 
Constant 
Constant 

Choice of perturbations including atmospheric 
retardation, perturbing planets, time-varying 
thrust orientation 

Thrust level & orientation, perturbation 
controls, position and velocity, ISP, a 
2-body initial conditions 
Not designed as optimization program 
Variable format 

High-thrust capture and departure, low-thrust 
spiralling 
Guidance analysis 

AVCO Know-How Document (internal documentation) 
Available for use on any funded study 
A 2-body version (1546D) is available. 


12 



Company/Organization 
Name of Program 
Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Specific Impulse 
Power Options 

Input (briefly) 

Output (briefly) 
Documentation 
Availability 
Other Comments 


Bell Telephone Laboratories 
DEPAK; MULTI 
A. J. Claus 
Company 

Description 

3 

3 

Cartesian 
Variable step 
Fortran IV 
IBM 7094, GE 635 
Variable, constant 

Constant, function of time, function of 
distance 

Position and velocity, Julian dates, ISP 
Time history 
Internal memoranda 

Written request addressed to Bell Labs required 
The program is essentially a general purpose 
numerical integration package with variable 
stepsize and operating in either single (8 
digits) or double (16 digits) precision. 
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Company/Organization 

Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Availability 
Other Comments 


: The Boeing Co. /Space Division - Guidance and 

Control 

: Variable Thrust Polynomial Optimizer 

: Forrester Johnson/Dari 1 Hahn 

: Company 

Description 

: 3 

: 2 

: Cartesian 

: Does not integrate 

: Fortran IV 

: Univac 1108 

: Flyby , rendezvous 

: Variable 

: Constant 

: Position and velocity, Julian dates, hyperbolic 

excess velocity, Isp 
: Self generating 

: /a 2 dt 

: Newton-Raphson 

: Time history, a 2 dt 

: Not included in program 

: Decided on an individual basis 

: Uses an approximate, but rapid, method 
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Comp any/Organizat ion 

Name of Program 
Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Other Options 

Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 

Other Applications 

Documentation 

Availability 


The Boeing Company/Space Div. , Guidance and 
Control Organization 

Space Trajectory Optimization Program (AS 2080) 

D. W. Hahn, Berdein F. Itzen 

Company 

Description 

3 

2 

Polar 
Fixed step 
Fortran IV 
Univac 1108 
Flyby, rendezvous 
Constant 

Constant, function of distance 
Planetocentric high-thrust departure to an 
input C 3 

Position and velocity, Isp, Vco 
Control schedule 

Coast times, thrust angle, final mass 

Gradient method 

Time history, payload 

High- thrust capture and departure; 

low- thrust spiralling 
General orbit transfer maneuvers, any 
thrust level 

Documented except for recent revisions 
Considered on an individual basis 
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Company/Organization 
Name of Program 
Authors/Contact 

Program Funded by 

Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 

Input (briefly) 


Initial Guesses Required 
Parameters Optimized 
Output (briefly) 

Planetocentric Phase 
Limitations of Program 
Documentation 
Availability 
Other Comments 


Reference 


Colorado State University 
Characteristic Length Method 
Professor W. R. Mickelsen, Dept, of 
Mech. Eng. 

NASA Grant NGR06-002-032 
Description 

Program does not compute trajectory. 

Fortran IV 
CDC 6400 

Flyby, rendezvous, round trip, etc. 

Variable, constant, other 
Constant, function of time, function of 
distance 

Hyperbolic excess speed and transfer time 
from impulsive thrust trajectories, launch 
vehicle payload versus V^, a of total 
electric propulsion system, Isp, T/W 0 
Propulsion time 

Isp, T/W 0 , Voo (by inspection of output) 
Approximate values of propulsion time, final 
mass of electric stage, final payload 
High- thrust capture and departure 
Trajectory similarity 
AIAA Paper No. 67-709 
Upon request 

The characteristic- length method is based 
on a nonrigorous assumption of similarity, 
between low-thrust and impulsive -thrust 
trajectories. Although the output is only 
approximate, it may serve well as a fast 
source for first guess inputs to complex 
programs . 

Zola, Charles L.: A Method of Approximating 

Propellant Requirements of Low-Thrust 
Trajectories. NASA TN D-3400, 1966. 
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Company/Organization 

Authors/Contact 


Number of Dimensions 
Number of Bodies 
Coordinate System 

Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Other Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Other Applications 
Other comments 


: Computer Applications, Inc. 

: Moshe Mangad, Ph.D. 

Description 

: 2, 3 

: 2, 3, 4 

: Cartesian, polar, nondimensional (i.e., Perkins 

Units) 

: Variable and constant 

: Fortran IV 

: IBM 7090/94, 704, Philco 5200 

: Flyby, rendezvous, moon transfer 

: Constant and variable 

: Constant, function of time, distance 

Coast period as f(t, d), f(t), or f(d) 

: Position and velocity, vehicle/payload 

constraints 
: T/W 

Isp, Voo, power level, T/W 0 
: Gradient method 

: Time history, Isp, T/W 

: Low thrust spiralling 

: Logistics study of Moon transfer 

4-body problem done for hyperbolic encounter 
study. 
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C omp any / Or g an i z at i on 

Name of Program 
Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Input (briefly) 
Parameters Optimized 
Other Applications 

Availability 

Other comments 


: University of Connecticut, Aerospace 

Engineering Dept . 

: PDER, PDERDP (double precision version) 

: Prof. Edward T. Pitkin 

: N. Conn., Research Foundation 

Description 

: 3 

: 2 

: Cartesian 

: N/A 

: Fortran IV 

: IBM 360/65 

: N/A 

Position and velocity at EPOCH 

: N/A 

: Computes transition matrix, inverse 2nd partials, 

perturbative derivatives, etc. 

: Single precision is share deck 3513; double 

precision version available from author. 

: Program is a general solution of two-body motion 

for all orbit shapes with 1st and 2nd partial 
derivatives plus perturbative derivatives for 
use in a variation of parameters integration. 
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Company/ Organization 
Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Other Options 

Initial Guesses Required 
Parameters Optimized 
(only p w in traj . 
part) 

Optimization Technique 
Output (briefly) 
Planetocentric Phase 

Documentation 

Availability 


General Dynamics/Convair 
Convair Low-Thrust Optimization Programs 
D. H. Higdon, J. F. Ingber, J. R. Grace 
Company 

Description 

2, 3 
2, N 

Cartesian 

Variable step, or fixed step 
Fortran IV 
CDC 6400 

Flyby, rendezvous, fixed time, free time 

Variable or constant 

Constant 

Complete optimization from launch to capture 
(excluding spirals) 

Lagrange multipliers 

Isp, T/W 0 , p w , Vo,, power level, coast time, 
thrust angle, final payload weight 

Newton-Raphson, maximum principle 
Time history, Isp, T/W, payload 
High-thrust capture and departure, low-thrust 
spiralling being developed 
Limited 
Proprietary 
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Company/Organization 
Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Documentation 
Availability 
Other Comments 


General Dynamics/Fort Worth Division 
Y28 Optimum Low-Thrust Procedure 
JPL/M. Poteet (Fort Worth Div.) 

Company 

Description 

2, 3 
2 

Polar 

Variable step 
Fortran IV 
IBM 7094 

Flyby, rendezvous 
Variable, constant 
Constant, function of time 
Position and velocity 
/a 2 dt 

Newton-Raphson 
/a 2 dt 

Not included in program 

JPL Engineer Planning Document #179 

Available from JPL 

JPL program is on production status at 
Fort Worth Division. 
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Company/Organization 
Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Other Comments 


General Dynamics/Fort Worth Division 
Low-Thrust Trajectory Optimization 
D . H . Krus e 
Company 

Description 

2, 3 
2 

Cartesian, polar 
Fixed step 
Fortran IV 
IBM 36- 
Rendezvous 
Constant 
Constant 

Julian dates, Isp, a 

Lagrange multipliers, flight time 

/a 2 dt, coast time, thrust angle 

Newton-Raphson 

/a 2 dt 

Not included in program 

Development of program is in initial phase 



Company/Organization 
Name of Program 

Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 

Parameters Optimized 
Output (briefly) 

Planetocentric Phase 
Limitations of Program 
Documentation 

Availability 


General Dynamics/Fort Worth Division 
Interplanetary Low-Acceleration Constant- 
Thrust Trajectory Patching Program 
S. W. Wilson/M. C. Poteet 
Company 

Description 

2 

2 

None 

Fortran IV 
IBM 7090 
Planetocentric 
Constant 
Constant 

Position, effective exhaust velocity, K p lanet > 
initial acceleration 

Initial (final) acceleration in heliocentric 
orbit 
None 

/a 2 dt, "effective time to escape" (or capture), 
AV, mass ratio 
Low- thrust spiralling 
10” 5 to 10 -2 local "g" acceleration 
Customer Utilization Instructions, Procedure 
V06 (Co. document) 

Proprietary, available under Government contract 
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Company/Organization 
Name of Program 

Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Documentation 

Availability 


General Dynamics/Fort Worth Division 
Low-Thrust Propulsion System Optimization 
Program 

S. W. Wilson/M. C. Poteet (Fort Worth Div.) 
Company 

Description 

2 

2 

None 

Fortran IV 
IBM 7090 

One way heliocentric 
Solved for in the program 
Constant 

Powerplant specific mass a, avg. acceleration 
/a 2 dt, efficiency coefficients 
None 

V '‘final- “initial- c (exhaust velocity), n 
Fibonacci an 

max V p w , y f . nal , a., c, n 
Not included in program 

Customer Utilization Instructions, Procedure 
V09 (Co. document) 

Proprietary; available under Government 
contract 
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Company/Organization 

Name of Program 
Authors /Contact 
Program Funded by 


Numbers of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 

Other Options 
Input (briefly] 
Parameters Optimized 
Optimization Technique 
Output (briefly] 
Planetocentric Phase 
Limitations of Program 
Documentation 
Availability 
Other Comments 


: General Electric Co. - Advanced Nuclear 

Systems Operation MSD 
: POLAR 

: Harold Brown 

: Company 

Description 

: 2 
: 2 

: Polar 

: No integration used 

: Fortran IV 

: IBM 7094 

: Flyby, rendezvous 

: Variable, constant 

: Constant, function of time, function of 

distance 
: See comment . 

: Position and velocity, Isp, a, Voo 

: T/W 0 , power level, coast time, thrust angle 

: Gradient method 

: Time history, payload, T/W 

: Not included in program 

: See comment . 

: Incomplete 

: Not recommended at this time 

: Program attempts to find a trajectory which will 

result in the specified power profile - 
constant or variable. Results to date have 
been unable to reduce deviations between the 
actual and specified power profiles to less 
than 10-15%. 
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C omp any / Or gan i z at i on 
Name of Program 
Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Other Options 

Input (briefly) 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 

Other Applications 
Documentation 
Availability 
Other Comments 


General Electric - Missile and Space Division 
Finite Bum Orbit Insertion Program 
D. Korenstein 
Company, NASA Contract 

Description 

3 

2 

Cartesian 
Fixed step 
Fortran IV 
IBM 7094 

Flyby, orbit insertion, escape 
Constant 

2-point boundary value problem for orbit 
insertion 

: Hyperbola, ellipse, engine, guidance, planet 

: T/W 0 

: Newton-Raphson 

Gravity loss 

: High-thrust capture and departure, low- thrust 

spiralling 

: Open-loop low-thrust trajectories 

: Not available 

Yes 

: Program is partially checked out (in process) 
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Comp any/Or gani zation 
Name of Program 
Contact 

Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 

Input (briefly) 

Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 


NASA Goddard Space Flight Center 
Variable Power, Variable Thrust 0ption/F369 
Kenneth I. Duck, Aux. Prop. Br. 

NAS 2-2928 

Description 

3 

2 

Cartesian 
Fixed step 
Fortran IV 
IBM 7094 

Flyby, rendezvous 
Variable 

Constant, function of distance, function of 
time 

Julian dates, hyperbolic excess velocity, 
powerplant specific mass, initial power 
/a 2 dt 

Newton- Raphson 

Time history, /a 2 dt 

Not included in program 
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Company/Organization 

Name of Program 
Contact 

Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 

Parameters Optimized 
Optimization Technique 
Output (briefly) 

Documentation 

Availability 


Goddard Space Flight Center, Mission and 
Trajectory Analysis Division 
Low-Thrust Trajectory Sub-Program 
Robert T. Groves 
NASA Contract NAS 5-10029 

Description 

2 or 3 
2 

Cartesian 
Fixed step 
Fortran IV 
IBM 7094 
Rendezvous 
Variable 
Constant 

Position and velocity, hyperbolic excess 
velocity 

Program uses ballistic trajectory for first 
fit 

/a 2 dt, thrust angle 
Newton -Raphson 

Time history of thrust acceleration and 
position, /a 2 dt 
In progress 

Program still under development 
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Comp any/Or gani zation 
Name of Program 

Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Other Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Other Applications 
D o cume nt at i on 
Availability 
Other Comments 


Grumman Aircraft, Bethpage, New York 
Low-Thrust Interplanetary Trajectory 
Optimization 

Hans K. Hinz and H. Gardner Moyer 
NASA-MSFC, Contract No. NAS 8-1549 

Description 

3 

2 

Cartesian 
Fixed step 
Fortran IV 
IBM 360/75 
Flyby, rendezvous 
Constant 

Thrust constant, thrust variable but limited 
Position and velocity and orbital elements 
of target planet 
Control function (s) 

Thrust angles and thrust magnitude 
Gradient method 
Time history 
Not included in program 
See Grumman report, RE -208, April 1965 
See final report of Contract NAS 8-1549 
Nonproprietary availability to interested users 
RE-208 discusses GEOCENTRIC MISSIONS AND USE OF 
MORE EFFICIENT GENERALIZED NEWTON -RAPHSON 
METHOD OF OPTIMIZATION. 
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Comp any /Organization 
Name of Program 
Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Other Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Documentation 
Other Comments 


Hayes International Corp. 

Orbital Elements for Low Thrust Trajectories 
Harry Passmore 
NASA Contract 

Description 

3 

2 

Orbital element 

Fixed step 

Fortran IV 

1620 in part 

Near Earth spiral 

Constant 

Constant thrust 

Position and velocity 

Lagrange multipliers 

Thrust angle 

Calculus of variations 

Time history 

Low -thrust spiralling 

NASA Report MTP-AERO-63-12 

Not presently operable 
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C omp any /Org an i z at i on 
Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Other Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Limitations of Program 
Documentation 
Availability 
Other Comments 


Hayes International Corp. 

Quasi-Steady Orbital Rendezvous 
Harry Passmore 
NASA Contract 

Description 

2 

Gravity of central body only 

Cartesian 

Fixed step 

Fortran IV 

IBM 360 

Transfer 

Constant 

Thrust magnitude constant 
Position and velocity, target position § 
velocity 
None 

Thrust angle 
Maximum principle 
Time history 
Not included in program 

Initial longitude between vehicle § target <90° 
Hayes Engineering Report 1324 
Computer program not published 
Only used to date for circular orbit 
rendezvous 
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Company/Organization 
Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 

Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 

Specific Impulse 
Power Options 
Other Options 


Input (briefly) 


Parameters Optimized 

Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Documentation 
Availability 


Hughes Aircraft Company 
Low-Thrust Trajectory 
D. MacPherson or J. H. Molitor 
Partially Supported Under JPL Contract 
No. 951144 

Description 

3 

N 

Output in radius, ecliptic longitude, and 
normal to ecliptic 
Variable step or fixed step 
Fortran IV 
GE 635 

Flyby; rendezvous; round trip can be simulated 
by inputting cases sequentially. 

Constant 

Function of time and/or function of distance 
Program converges automatically to a trajec- 
tory which satisfies terminal conditions; 
optimization is by parameter search. 

Isp, powerplant specific mass, thrust 
profile, spacecraft mass, power profile, 

Vo,, launch date 

Program converges automatically to a trajec- 
tory which satisfies terminal conditions. 
Parameter search 
Time history 

High- thrust capture, and departure 

Not complete 

Subject to negotiation 
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Company/Organization 
Name of Program 
Authors /Contact 

Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 

Other Applications 

Documentation 

Availability 


IBM Federal Systems Division 
OPGUIDE 

G. W. Johnson, 1730 Cambridge Street 
Cambridge, Mass. 

NASA Contract 

Description 

3 

2 

Cartesian 
Variable step 
Fortran IV 

IBM 7094 - IBM System 360 

Rendezvous 

Variable 

Function of time 
Position and velocity, Isp 
Lagrange multipliers 
Coast time, thrust angle 
Newton-Raphson 
Time history 

High thrust orbital transfers 

Preliminary 

Limited 
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Company/Organization 
Name of Program 
Author s /C ont act 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 

Parameters Optimized 
Optimization Technique 
Output (briefly) 

Other Applications 
Documentation 


IIT Research Institute 

LTNAV-1 (Low-Thrust Navigation Analysis) 

Alan L. Friedlander 

NASA Contract 2-22401 and 2-3348 

Description 

3 

2 

Cartesian, polar 
Variable step, fixed step 
Fortran II 
IBM 7094 

Flyby, rendezvous 
Variable, constant 
Constant 

Position and velocity 

Lagrange multipliers for converged optimal 
trajectory 
Thrust angle 
Calculus of variations 
Time history 
Navigation studies 
NASA CR-457 
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Company/Organization 
Name of Program 

Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Other Options 
Input (briefly) 

Output (briefly) 
Planetocentric Phase 
Limitations of Program 
Documentation 
Availability 


: I IT Research Institute 

: LTNAV-2 (Low-Thrust Planetocentric 

Navigation Analysis) 

: ALAN L. Friedlander 

NASA Contract 2-2401 and 2-3348 

Description 

: 3 

: 3 

: Cartesian 

: Variable step, fixed step 

: Fortran IV 

: IBM 7094 

: Planetocentric only 

Constant 
: Constant 

: Constant acceleration 

: Position and velocity, orbit elements 

: Time history 

Program treats planetocentric phase only 
: Tangential thrust only 

: NASA CR 740 

: Yes, but documentation is negligible. 
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Company/Organization 
Name of Program 

Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 

Specific Impulse 
Power Options 
Input (briefly) 


Initial Guesses Required 
Parameters Optimized 

Optimization Technique 
Output (briefly) 
Planetocentric Phase 


Other Applications 
Documentation 

Availability 


: JPL 

Optimum Thrust Programs for Power Limited 
Propulsion Systems (1016000) 

: Carl Sauer 

: NASA Contract 

Description 

: 2, 3 

: 2 

: Polar 

Variable step, fixed step 
: FAP 

: IBM 7094 

: Flyby, rendezvous, solar probe, out of 

ecliptic, maximum energy 
: Variable, constant 

: Constant, function of distance 

: Position and velocity, Julian dates, 

hyperbolic excess velocity, Isp, power- 
plant specific mass, a 
: Lagrange multipliers, T/W 

: Isp, T/W 0 , p w , Vco, power level, fa 2 dt, 

coast time, thrust angle 
: Newton-Raphson 

: Time history, Isp, T/W, payload, /a 2 dt 

High-thrust capture and departure, low-thrust 
spiralling; approximations made to helio- 
centric portion using "asymptotic velocity 
bias." Spirals and high thrust escape can 
be run separately. 

: Vehicle path optimization in vacuum 

: JPL EPD 179 11/1/63, latest documentation 

lacking 

: Earlier version being used 



C omp any/Organization 
Name of Program 
Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 


Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Limitations of Program 
Initial Guesses Required 
Availability 


JPL 

ASTRAL 

D. Alderson/W. Stavro 
NASA Contract 

Description 

3 (limited 3rd dimension) 

2 

Spherical 

Analytic approximation 
Fortran IV 
IBM 7094 

Flyby, rendezvous 
Constant 

Function of distance 

Julian dates, Isp, powerplant specific mass, 
planet (target, launch), powerplant 
efficiency 

Isp, T/W q , y , Vo,, coast time, thrust angle 

"Complex” method 

Time history, Isp, T/W, payload 

High -thrust capture and departure 

Travel angle, T/local g, eccentricity of orbit 

None 

Not at present 
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Company/Organi zation 
Name of Program 
Auth or s / C ont act 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 


Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Documentation 


JPL 

MALTS 

TRW Systems/J. Driver 
NASA Contract NAS 7-503 

Description 

3 

Up to 10 
Spherical 

Cowell or Adams Moulton 
Fortran IV 
IBM 7094 

Flyby, rendezvous 
Constant 

Function of distance 

Position and velocity, Julian dates, excess 
velocity, Isp, powerplant specific mass, 
thrustor efficiency, time of flight 
Lagrange multipliers, or orbit plane 
fa 2 dt, time, final mass, payload 
Calculus of variations 
Time history, /a 2 dt 
User's Manual available 
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Company/Organi zat ion 
Name of Program 
Authors/Contact 
Program Funded by 


NASA Langley Research Center 
PRESTO (DELTA) 

Robert E. Smith 
NASA Contract 

Description 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Other Options 

Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Other Applications 
Documentation 
Availability 
Other Comments 


3 

2 

Polar 
Fixed step 
Fortran IV 
CDC 6600 
Rendezvous 
Constant 
Constant 

Earth to orbit transfer. Earth to Lunar 
transfer 

Position and velocity 

Control history, coast angle 

Coast time, thrust angle 

Gradient method 

Time history 

Not included in program 

Earth orbit missions. Lunar missions 

NASA CR-158 


See note below • 

Program applicable for preliminary studies 
NOTE: The Presto Program is available from the COSMIC Computer Center, 

University of Georgia, Athens, Georgia, 30601. Their identification 
number for the program is LAR- 10078. 


38 



Comp any /Organization 
Name of Program 

Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Documentation 
Availability 
Other Comments 


NASA Langley Research Center 
TOMSD (Trajectory Optimization by Method of 
Steepest Descent) 

Robert E. Smith 

Air Force-NASA contract 

Description 

3 

1 or more 
Rectangular 
Variable or fixed step 
Fortran IV 
CDC 6000 

Reentry, boost, rendezvous 
Position and velocity 
Control history 

Angle of attack, bank angle, initial state 

Gradient method 

Time history 

Not included in program 

Technial Report AFFDL-TR-67-108, Vols . I, II, III 
Upon request 

This program is applicable to a variety of 
trajectory optimization problems. 
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C omp any /Organization 
Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Other Options 
Input (briefly) 

Initial Guesses Required 

Parameters Optimized 

Optimization Technique 
Output (briefly) 
Planetocentric Phase 

Other Applications 
Documentation 
Availability 
Other Comments 


NASA Lewis Research Center 
Lewis N-Body Code 
William C. Strack 
Company 

Description 

2, 3 
2, 3, 8 

Cartesian, orbit elements 
Variable step, fixed step 
Fortran IV 

IBM 7094, CDC 6600, IBM 7044 
Flyby, rendezvous 
Constant 

Constant, function of distance 
Targeting iteration, choice of thrust control 
Position and velocity, Julian dates, Isp, a, V ro , 
tank and structure fraction 
Lagrange multipliers, T/W, thrust angle 
n, 4, etc. 

Isp, T/W 0 , Voo, power level, coast time, 
thrust angle 

Calculus of variations, univariate search 
Time history, payload 

High-thrust capture and departure, low-thrust 
spiralling 

Booster ascent paths 

TN D-1730 with supplementary material 

Yes: 25 copies already sent out. 

Several versions of code exist (see 
NASA TM X52322) . 


40 



LOW-THRUST TRAJECTORY AND MASS COMPUTATION PROGRAM 


Company Organization 
Name of Program 
Authors /Contact 
Program Funded by 


NASA-Lewis Research Center 
MAXIM 

E. A. Willis, Jr. 

Lewis Research Center (in-house) 
Description 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 

Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Other Options 
Input 

Initial Guesses Required 
Parameters Optimized 


Optimization Technique 
Output 

Documentation of Program 

Availability to Other 
Users 


2 

2 

Polar 

Runge-Kutta with variable step size and error 
control 
Fortran IV 
IBM 7094 

Orbit escape, capture, or transfer 
Constant I 

Constant power implied 
Constant thrust, constant acceleration 
Position and velocity; vehicle design parameters 
for each stage (up to 5) 

Initial positions and one multiplier (X Q ) 

Thrust steering and switching program, staging 
points and acceleration levels are optimized 
for minimum AV or weight 
Maximum principle with conjugate-gradient search 
for the 2-point b.v.p.'s 
Time history, AV's, weights 
Theory and typical application covered in TN 
D-3606 , TN D-4534 , and TN D-5011. 

Considered individually. Operating instructions 
are in the form of internal notes only. 
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Comp any /Organi zation 

Name of Program 
Authors / Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 


Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 

Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Limitations of Program 
Documentation 


Availability 


: Research Laboratories; Lockheed Missiles 

and Space Company 

: Low-Thrust Heliocentric Transfers 

: H . E . Rauch 

: Company 

Description 
: 2 or 3 

: Vehicle plus up to three planets 

: Cartesian 

: Fixed step 

: Fortran IV 
: Univac 1108 

: Flyby, rendezvous, and round trip with 

variable acceleration only, rendezvous with 
constant acceleration 
: Variable, constant 

: Constant 

: Position and velocity of planets, elapsed 

time 

: Lagrange multipliers 

: /a 2 dt, coast time, thrust angle for constant 

acceleration 

: Guessing adjoint variables 

: Time history, /a 2 dt 

: Not included in program 

: Not "production" type program 

: "Low Thrust Swingby Trajectories," by 

H. E. Rauch, XVIIIth International Astro- 
nautical Congress - Belgrade, 1967, Pergamon 
Press, Oxford, New York, PWN-Polish Scientific 
Publishers, 1968, pp. 269-284. 

: Could be made available 
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Company/Organization 

Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Specific Impulse 
Power Options 

Input (briefly) 


Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Limitations of Program 
Other Applications 

Documentation 


Availability 


Research Laboratories, Lockheed Missiles 
and Space Co. 

Low-Thrust Escape From Planet 

H . E . Rauch 

Company 

Description 

2 

2 

Both, but primarily polar 
Variable step 
Fortran IV 
Univac 1108 

Variable, constant acceleration 
None, trajectory with constant acceleration 
in velocity direction 

Initial position and velocity, final desired 
hyperbolic excess velocity, constant or 
"average" acceleration 
Lagrange multipliers 

fa 2 dt, thrust angle for constant acceleration 
Guessing adjoint variables 
Time history, fa 2 dt 
Low-Thrust spiralling 
Not "production" type program 
Uses asymptotic matching to calculate 
heliocentric "offsets" 

"Optimum Guidance for a Low Thrust Interplane- 
tary Vehicle," J. Breakwell and H. Rauch, 
AIAA Journal, Vol. 4, No. 4, April 1966 
Could be made available 
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Company/Organization 

Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 

Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Limitations of Program 
Documentation 

Availability 
Other Comments 


Lockheed Missiles § Space Co. /Advanced 
Flight Mechanics 
LTC0MP 

M. A. Krop, Dept. 55-34, B. 201, Ext. 45869 
Company 

Description 

2, 3 
4 

Planet-center § aligned along the excess 
velocity vectors 
None 

Fortran IV 
Univac 1108 

Single- leg interplanetary transfer 

Variable, constant 

Constant 

Julian dates, planet numbers 
Time offsets at terminal planets 
/a 2 dt 

Analytic approximation 
/a 2 dt 

Not included in program 
Travel angle 

Method described in Chapter 4 of LMSC 
Rept. No. 3-17-64-1, Apr. 30, 1964 
Still being developed 

LTC0MP rapidly simulates the propulsion 
requirements for a low-thrust transfer. 
Although quite simplified, it computes 
optimum /a 2 dt to within a couple of 
percent of rigorously optimized values. 
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Company/Organization 
Name of Program 


Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 


Los Alamos Scientific Laboratory, Los Alamos, 
New Mexico 

1965 Version, produced by Dept, of Aerospace 
and Mechanical Sciences, Princeton Univ., 
Report No. 7170-1 
AEC 

Description 

2 

2 

Polar 

Variable step 

Fortran IV 

IBM 7094, IBM 7030 

Rendezvous 

Constant 

Constant (with coast) 

Position and velocity, Isp 
Lagrange multipliers 
Coast time, thrust angle 
Direct 

Time history, T/W 
Not included in program 
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Company/Organization 
Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Input (briefly) 

Parameters Optimized 
Optimization Technique 

Output (briefly) 
Planetocentric Phase 
Limitations of Program 


McDonnell Douglas Astronautics Company 

FYNDIF 

F . I . Mann 

Company 

Description 

2 

2 

Cartesian 

Second order finite difference procedure 
Fortran IV 
Univac 1108 
Rendezvous, intercept 

Position and velocity, hyperbolic excess 
velocity 
: / a 2 dt 

: Calculus of variations and Van Dine's 

Generalized Newton-Raphson 
: Time history 

: Not included in program 

: Experimental program, planets in circular 

orbits, few options 


46 



Comp any/Organization 
Name of Program 

Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 

Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Other Applications 
Documentation 
Availability 


McDonnell Douglas Astronautics Company 
LOWTOP, High-Low Thrust Trajectory Optimization 
Program 

Jerry L. Horsewood/F. I. Mann 
Company 

Description 

2 

2 

Cartesian 

Variable step, fourth-order Runge Kutta 

Fortran IV 

Univac 1108 

Flyby, rendezvous 

Variable, constant 

Constant or variable 

Position and velocity, Julian dates, V m> Isp, 
thrust, power, weight-flow 
Lagrange multipliers 

/a 2 dt, coast time, thrust angle, flight time, 
final speed 

Newton-Raphson, maximum principle 
Time history, Isp, T/W, velocity loss 
High -thrust capture and departure 
See company documentation on LOWTOP 
First draft complete 
Currently not available 
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Company/Organization 
Name of Program 
Authors /Contact 
Program funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 

Other Applications 


Documentation 
Availability 
Other Comments 


McDonnell Douglas Astronautics Company 
WHAMO, Weight History and Mission Optimizer 
F . I . Mann 
Company 

Description 

2 

2 

Cartesian 

Variable step, Fourth-order Runge Kutta 
Fortran IV 
Univac 1108 

Rendezvous, round trip, multi-leg missions 
Variable or constant 
Variable or constant 

Julian dates, hyperbolic excess velocity, Isp, 
thrust, power, weight flow 
Lagrange multipliers, T/W 

Isp, coast time, thrust angle, gross weight 
Newton-Raphson 
Time history, Isp, T/W 
High capture and departure with engine 
shutdown § restart 

Simulates both high and low thrust engine 
shutdown and restart with high thrust 
staging according to the non-integral-burn 
method 
Partial 

Currently not available 
Program not extensively used yet 
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Comp any /Organization 
Name of Program 
Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Input (briefly) 

Initial Guesses Required: 
Parameters Optimized 

Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Limitations of Program 
Other Applications 

Documentation 

Availability 


McDonnell Douglas Astronautics Company 
HITOP 

F . I . Mann 
Company 

Description 

3 

2 

Cartesian 

Impulse trajectories analytic solution 

Fortran IV 

Univac 1108 

Rendezvous 

Julian dates 

Time and place of deep space impulse 
Either (V«>i + AV 2 + V 003 ) or 
(AVi + AV2 + AV3) 

Newton-Raphson or transversality procedure 
Primer vector, V™, AV, orbit parameters 
High -thrust capture and departure 
Currently only one deep-space impulse 
Program may be used to study the high-thrust 
limit of low-thrust trajectories. 

: None 

: Currently not available 



Company/Organi zation 
Name of Program 

Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language: 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Limitations of Program 
Other Applications 
Documentation 
Other Comments 


: Space Division, North American Rockwell Corp. 

Optimum Trajectory Program for Power- Limited 
Propulsion Systems 
: Russell P. Nagorski 

: Company 


Description 


3 

4 

Polar 
Fixed step 
I II 

IBM 7094 

Flyby, rendezvous 
Variable, constant 
Constant 

Position and velocity, Julian dates, V , Isp, a 
Lagrange multipliers, T/W 
/a 2 dt, coast time 
Newton-Raphson 

Time history, Isp, T/W, payload, /a 2 dt 
Not' included in program 
Initial guess, iteration logic 
General intercept and rendezvous problems 
User's Manual complete 

Derivative of JPL (Melbourne and Sauer) program 
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Company/Organization 
Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 

Documentation 


: Mission Analysis Division, OART/NASA 

: BABE 

: A. C. Mascy/S. W. Pitts 

: In-house 

Description 

: 2 
: 2 
Polar 

: Variable step and fixed step 

Fortran IV 
: IBM 7094 

: Flyby, rendezvous 

: Constant 

: Constant 

: Julian dates, hyperbolic excess velocity, a 

: Isp, y^, coast time 

Gradient method 
Coast time, Isp, payload, y^ 

High-thrust capture and departure, low-thrust 
spiralling 
: In progress 
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Company/Organization 
Name of Program 
Authors /Contact 
Program Funded by 


Mission Analysis Division, OART/NASA 
QUICKLY 

Alfred C. Mascy 
In-house 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 


Program Language 
Computer Now Used 
Trajectory Options 

Specific Impulse 
Power Options 
Other Options 


Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 

Planetocentric Phase 

Other Applications 

Documentation 

Availability 


Description 

Program uses functional relationships for the 
energy requirements of pre-computed built-in 
performance data which were obtained by accu- 
\ rate 3-dimensional, 2-body, Cartesian 

coordinate, Newton-Raphson finite-difference 
programs . 

: Fortran IV 

: IBM 360/50 

: Flyby, rendezvous, solar probe, extra-ecliptic, 

comet intercept 
: Constant 

: Constant, function of time or distance 

: Optimum or fixed power level or alpha as a 

function of power. Also, optimum or fixed 
hyperbolic excess velocities 
: Trip time, power plant specific mass a, launch 

vehicle name, capture and departure mode 
: None 

: Isp, y^, Vcoi, Vcc> 2 , power level, coast time 

: Gradient method and Newton-Raphson 

: Payload, Isp, y^, Voo 1} Voo 2 , power level, powered 

time, also, plot option 

: High-thrust capture and departure, low-thrust 

spiralling, launch vehicle boost 
: May be used to determine performance of two- 

impulse ballistic systems for comparisons 
: In progress 

: Upon request 
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C omp any /Or g an i z at i on 
Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 

Other Options 
Input (briefly) 


Initial Guesses Required 
Parameters Optimized 

Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Other Applications 

Documentation 
Availability 
Other Comments 


ASMAR, Princeton University 
Lion 1 

J. Campbell/P. M. Lion 
NASA contract 

Description 

3 

2 

Cartesian 
Variable step 
Fortran IV 

IBM 7094, CDC 6600, IBM 360/50 

Flyby, rendezvous, round trip, escape* 

Constant 

Constant, function of time, function of 
distance 

Booster selection 

Position and velocity or Julian dates, power 
plant specific mass, tankage, structure 
factors 

Lagrange multipliers 

Isp, T/W , Voo, power level, coast time, 
thrust angle* 

Neighboring extremal 
Time history, Isp, T/W, payload 
High- thrust capture and departure 
Iterator can be used to solve any nonlinear 
(algebraic) equations. 

In progress 
Available on request 

*0ther options available depending on input 
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Company/Oxganization 
Name of Program 
Authors /Contact 
Program Funded by 


Number cf Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 

Other Options 
Input (briefly) 

Initial Guesses Required 
Optimization Technique 
Planetocentric Phase 
Other Applications 
Documentation 
Availability 
Other Comments 


ASMAR, Princeton University 
Campbell 1 (Modified Item) 

J. Campbell/P. M. Lion 
NASA contract 

Description 

3 

9 

Cartesian 
Variable step 
Fortran IV 

IBM 7094, IBM 360/50 
Flyby, rendezvous 
Constant 

Constant, function of time, function of 
distance 

Programmed thrust 

Position and velocity or Julian dates 

Lagrange multipliers 

Neighboring extremal 

Included in program 

Uses generalized iterator 

Partially complete 

Difficult to transfer 

Used mainly as check on 2-body program 
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Comp any /Organization 
Name of Program 
Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 

Trajectory Options 

Specific Impulse 
Power Options 

Other Options 
Input (briefly) 


Initial Guesses Required 
Parameters Optimized 


Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Documentation 
Availability 
Other Comments 


ASMAR, Princeton University 
Gordon 2 

C. N. Gordon/P. M. Lion 
NASA contract 

Description 

2 

2 

Polar 

Fixed step 
Fortran IV 

IBM 7094, CDC 6600 (now being adapted), 

IBM 360/50 

Flyby, rendezvous, round trip, constrained 
velocity flyby 
Constant 

Constant, function of time, function of 
distance 
Select booster 

Position and velocity or Julian dates (can 
input or optimize: excess velocity, Isp) , 

power plant specific mass, a, tankage, 
structure factors 
Lagrange multipliers 

Isp, T/W 0 , y w , Voo, power level, direction 
of arrival velocity, coast time, thrust 
angle, C 3 

Neighboring extremal 

Time history, Isp, T/W, payload, summary 
High -thrust capture and departure 
In progress 
Available on request 

Also being used at Battelle, IITRI, MSFC 
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Company/Organization 
Name of Program 
Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Input (briefly) 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Documentation 
Availability 
Other Comments 


ASMAR, Princeton University 
PRIMER 

M. Minkoff/P. M. Lion 
NASA contract 

Description 

3 

2 

Cartesian 

Impulsive trajectories analytic solution 
Fortran IV 
IBM 7094 

Flyby, rendezvous 

Position and velocity or Julian dates 

Conjugate gradient 

Time history, J = 2 ^ ] A V ] 

High- thrust capture and departure 
In thesis form 
Available on request 

Finds optimum impulsive trajectory regardless 
of the number of impulses 
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Company/Organization 
Name of Program 
Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Other Options 

Input (briefly) 

Output (briefly) 

Planetocentric Phase 

Limitations of Program 

Other Applications 
Documentation 

Availability 


: Rand Corp . 

: ROCKET 

: B. W. Boehm 

: USAF Project RAND 

Description 

: 3 

: 4 

: Cartesian, polar 

: Variable step, fixed step, Cowell method 

: Fortran IV, II 

: IBM 7094, CDC 6600, Univac 1108, GE 635 

Flyby, departure and capture 
: Variable, constant 

: Constant, function of time 

: Numerous propulsion, aerodynamic, guidance 

options 

: Position and velocity 

Time history, Isp, T/W, osculating orbit 
parameters, tracking obs . 

: High- thrust capture and departure, low-thrust 

spiralling 

: Cowell method makes ROCKET inefficient for 

long-duration trajectories. 

: Boost, reentry, airbreathing propulsion 

: Book: ROCKET, by B. W. Boehm, Prentice-Hall, 

1964 

: Through .SHARE Distribution Agency, 112 E. Post 

Road, White Plains, New York (send blank mag 
tape, ask for #3001 (Fortran II version) 
or #3485 (Fortran IV version)) 



C omp any /Or gan i z at i on 
Name of Program 
Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 

Other Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 

Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Documentation 
Availability 


RCA/Astro Electronics Division 
Gordon 2 

C. N. Gordon/ J. M. L. Holman 
Company, NASA contract 

Description 

2 

2 

Polar 

Fixed step 
Fortran IV 

IBM 7094, IBM 360/50 
Flyby, rendezvous 
Constant 

Constant, function of time, function of 
distance 

Choice of launch vehicle, sweep of parameters 
Julian dates, power plant specific mass 
Lagrange multipliers, T/W, power, C 3 
Isp, T/W Q , power level, coast time, 
thrust angle, C 3 
Neighboring extremal 
Time history, Isp, T/W, payload 
High-thrust capture and departure 
In progress 

Available through ASMAR, Princeton University 
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Corap any/ Organization 

Name of Program 

Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Other Applications 

Documentation 
Availability 
Other Comments 


System Sciences Corp. (Subsidiary of Computer 
Science Corp.) 

High-Precision, Low-Thrust Integration and 
Optimization Prog. 

Pitkin and Baker/Dr. R. M. L. Baker, Jr. 
Company 

Description 

3 

Any number for which ephemerides are given 

Cartesian 

Variable step 

II 

IBM 7094 
Open 

Variable 
Function of time 
Position and velocity, Isp 
Lagrange multipliers 
Power level 
Gradient method 
Time history 
Low-thrust spiralling 

Can interface with other impulsive thrust 
programs 

Listing of basic program 
Yes 

Utilizes universal variable so that it can 
pass evenly from ellipses through parabolas 
to hyperbolas 



Company/Organization 
Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 

Limitations of Program 

Documentation 

Availability 


U. of Texas, Austin, Dept, of Aerospace Eng. 
NMS 

G. J. Lastman § W. T. Fowler 

NASA contract and graduate student and faculty 
time contributions 

Description 

3 

2 

Cartesian 
Fixed step 
Fortran IV 
CDC 6600 
Rendezvous 
Constant 
Constant 

Position and velocity, or Julian dates 
Lagrange multipliers 
/a 2 dt, flight time 
Perturbation 

Time history, initial multipliers 
Not inluded in program, capture being 
implemented 

Travel time limit of 1000 days 
Limited 

Listing and/or deck on request 
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Comp any/Organizat ion 
Name of Program 
Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Porgram Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 

Other Options 

Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 

Limitations of Program 
Documentation 

Availability 
Other Comments 


TRW Systems/Analytic Mechanics Dept. 

MALTS, Mission Analysis § Simulation Program 
D. B. Smith fj W, D. Dickerson (R3/T, 3442) 
NASA contract 

Description 

3 

N 

Cartesian 

Variable step, fixed step (both) 

Fortran IV 
IBM 7094 

Flyby, rendezvous 
Variable, constant 

Constant, function of time, function of 
distance 

Nonoptimal or optimal steering, computed or 
tabular thrust, others, see report 
Position and velocity or Julian dates, V m , 
Isp, a 

Steering angles 

Can optimize any selected parameter 
Newton-Raphson (modified) 

Time history, Isp, T/W, payload 
High-thrust capture and departure 
Optimal phasing logic not included yet 
User's Guide (Vol. I § II), Interim Report 
NAS 7-503 

See John Driver, JPL 

The program is not completely ready for 
general use. 
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Company/Organization 
Name of Program 

Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Limitations of Program 
Documentation 
Availability 
Other Comments 


United Aircraft Research Laboratories 
Mass ratio optimization (simplified § improved 
payload fraction definition) F487 
R. Ragsac/Betty Knose 
Company 

Description 

3 

2 

Cartesian 
Fixed step 
Fortran V 
Univac 1108 
Rendezvous 
Constant 
Constant 

Julian dates, a, n, thrustor specific mass 
curve 
Isp, ui 

Isp, y w , /a 2 dt, coast time 
Newton-Raphson 

Isp, /a 2 dt, y , yi, y T powered time 
Not included in program 

Travel angle, powered time is approximate 

Sample input and explanation 

Yes 

Differs from F530 in capability to maximize 
realistic definition for payload fraction 
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Comp any/Organization 
Name of Program 
Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Documentation 
Availability 
Other Comments 


United Aircraft Research Laboratories 
High § Low Thrust Mass 0ptimization/F211 
Betty Knose/R. Ragsac 
NASA Contract NAS 2-2928 

Description 

2 or 3 
2 

See other comments 
See other comments 
Fortran V 
Univac 1108 

Flyby, rendezvous, round trip 
Constant 

a, high and low-thrust propulsion parameters, 

table of y T , p Isp vs. 

L w c 

None 

Mass on Earth parking orbit 
Search 

System masses and propulsion parameters 
Hyperbolic excess speeds (high-thrust system) 
User's Manual 
See NASA, MAD 

Mass optimization only; trajectory data 
are input 



Company/Organization 
Name of Program 

Authors/Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 

Planeto centric Phase 
Limitations of Program 
Documentation 
Availability 
Other Comments 


United Aircraft Research Laboratories 
Mass Ratio Optimization (simplified payload 
fraction) F530 
R. Ragsac/Betty Knose 
Company 

Description 

3 

2 

Cartesian 
Fixed step 
Fortran V 
Univac 1108 
Rendezvous 
Constant 
Constant 

Julian dates, a 
None 

Isp, u w , /a 2 dt, coast time 
Newton-Raphson 

Isp, /a 2 dt, y , y T , pi powered time 
Not included in program 

Travel angle ~3tt, powered time is approximate 

Sample input and explanation 

Yes 

Optimization of Isp and y^ uses accurate 
approximation technique, single coast 
period only 
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Company/Organization 
Name of Program 
Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 

Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 
Planetocentric Phase 
Documentation 
Availability 


United Aircraft Research Laboratories 
Variable Power, Variable Thrust 0ption/F369 
Mrs. E. Knose, Mr. G. Thrasher, UARL 
NASA Contract NAS 8-11309, NAS 2-2928 

Description 

3 

2 

Cartesian 
Fixed step 
Fortran IV 

Univac 1108, Fortran V 

Flyby, rendezvous, solar probe 

Variable 

Constant, function of time, function of 
distance 

Position and velocity, Julian dates, Voo, power 
option input 
None 
/a 2 dt 

Newton-Raphson 
Time history, /a 2 dt 
Not included in program 
Sample input and explanation 
Determined by NASA 
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Company/Organization 
Name of Program 
Authors /Contact 
Program Funded by 


Number of Dimensions 
Number of Bodies 
Coordinate System 
Integration Technique 
Program Language 
Computer Now Used 
Trajectory Options 
Specific Impulse 
Power Options 
Other Options 

Input (briefly) 

Initial Guesses Required 
Parameters Optimized 
Optimization Technique 
Output (briefly) 

Planetocentric Phase 
Limitations of Program 
Documentation 
Availability 
Other Comments 


United Aircraft Research Laboratories 
Constant-Thrust, Multiple-Coast/F615 
C. P. Van Dine/Betty Knose 
NASA Contract NAS 2-2928 

Description 

2 or 3 
2 

Cartesian 
Fixed step 
Fortran V 
Univac 1108 
Flyby, rendezvous 
Constant 
Constant 

Optimize Isp and y^, fix y^ opt Isp, fix Isp 
opt y , or input Isp and y^ 

Julian dates, V^, a, efficiency parameter 
None 

Isp y , /a 2 dt, coast time 
r w 

Newton-Raphson 

Time history, Isp, payload, fa 2 dt, powered 
times 

Not included in program 
Travel angle ~3 tt 

Programmer's Manual, UARL F-910352-13 
See NASA, MAD 

Two heliocentric coast periods allowed 


National Aeronautics and Space Administration 
Moffett Field, Calif., 94035, May 17, 1968 
130-0604-07 
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